#include <iostream>
#include "Matrix4.h"
#include "Vector.h"
using namespace std;

float angulo_entre_vetores(Vector3 v, Vector3 u)
{
    return(acos((v^u)/(v.length()*u.length())) * (180/M_PI));
}

int main(void)
{
    Point3 A(1,2,0);
    Point3 B;

    Matrix4 SH = shearMatrix(45);

    cout<<"SH:"<<endl;SH.imprimir();
    cout<<"\nA = ";A.imprimir();
    B = SH * A;
    cout<<"\nA' = ";B.imprimir();

    cout<<"\nAngulo A'-(Ax,0,0) e eixo-X = "<<angulo_entre_vetores(B-Vector3(A(0),0,0),Vector3(1,0,0))<<endl;

    return(0);
}
